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1 .  INTRODUCTION 


The  muzzle  blast  associated  with  the  firing  of  a gun  can  be  of 
sufficient  magnitude  to  cause  considerable  damage  to  structures  or 
personnel  some  distance  away  from  the  muzzle.  Analytical  methods  of 
predicting  the  characteristics  of  the  blast  field  are  of  primary  im- 
portance to  the  analysis  of  a weapon  system  to  insure  structural 
integrity  and  personnel  safety.  Ideally,  the  analytical  methods 
should  be  based  on  basic  principles  so  that  the  blast  field  can  be 
calculated  without  resort  to  empirically  determined  parameters. 

Westine"^  has  applied  the  technique  of  dimensional  analysis  to 
the  blast  field  and,  based  on  the  examination  of  large  quantities  of 
experimental  data,  he  presents  empirical  formulas  for  the  prediction  of 
blast  overpressure,  pulse  length  and  time-of-arrival . The  numerical 
constants  in  his  equations  are  given  for  several  barrel  elevation 
angles  together  with  a coefficient  which  determines  the  effective 
energy  release  for  several  propellents. 

A completely  different  approach  for  the  prediction  of  blast 
fields  has  recently  been  devised  by  Henriksen  and  Cummings^.  Their 
analysis  utilizes  the  Unified  Theory  of  Explosions  (UTE)  developed  by 
F.  B.  Porzel^  in  the  course  of  an  extensive  study  of  nuclear  and  HE 
explosions.  In  their  theory  Henriksen  and  Cummings  (HC)  apply  basic 
principles  to  evaluate  the  fraction  of  propellent  energy  which  is 
available  to  the  blast  field  and  then  use  UTE  to  describe  the  propa- 
gation of  the  blast  with  distance.  A function  was  derived  to  correlate 
the  nonspherical  muzzle  blast  field  to  the  spherical  blast  field 
predictions  of  UTE. 

In  the  present  study  the  basic  approach  of  HC  has  been  retained, 
but  the  method  of  calculating  the  equivalent  hydrodynamics  yield  has 
been  modified;  the  non-spherical  nature  of  the  gun  muzzle  blast  field 
has  been  included  in  UTE  and  a completely  new  method  of  calculating 
pulse  length  has  been  developed. 


1.  P.  S.  Westine  and  J.  C.  Hokanson,  "Prediction  of  Stand-Off  Distances 

to  Present  Loss  of  Hearing  from  Muzzle  Blast" , R-CR-7 5-003  South- 
west Research  Institute,  February,  1975.  Ad/A-005  274. 

2.  B.  B.  Henriksen  and  B.  E.  Cummings,  "An  a priori  Theory  for  Muzzle 
Blast  Overpressure  and  Pulse  Length  Determination" , BRL 

3.  F.  B.  Porzel,  "Introduction  to  a Unified  Theory  of  Explosions (UTE) " , 
NOL  72-209 , US  Naval  Ordance  Laboratory , White  Oak,  Silver  Spring, 
MD,  September,  1972.  AD  758000. 
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II.  THEORY 

A.  Porzel’s  Unified  Theory  of  Explosions 

Both  HC  and  the  present  development  depend  heavily  on  the  theory 
and  computational  methods  of  Porzel's  UTE-^.  This  theory  was  developi-.i 
over  a period  of  time  and  appeared  in  a number  of  company  and  lab- 
oratory reports.  A fairly  comprehensive  review  of  this  theory  appears 
in  Reference  2 and  only  a brief  synopsis  will  be  presented  here. 

The  UTE  Is  strongly  dependent  on  the  separation  of  the  total 
available  energy  into  "prompt  energy"  which  is  avail-.ble  to  drive  the 
blast  and  delayed  energy.  The  prompt  energy  includes  the  energy  due 
to  static  overpressure  and  dynamic  pressure  and  also  includes  the 
kinetic  energy  due  to  material  velocity.  The  remaining  energy  in  the 
explosive  is  delayed,  meaning  that  it  is  transported  too  slowly  to 
support  the  blast. 

This  division  of  energy  removes  some  of  the  mathematical  diffi- 
culties associated  with  the  initial  instant  of  energy  release  and 
permits  an  accurate  calculation  of  the  prompt  energy  remaining  in  the 
blast  wave  as  it  spreads  outward  and  decreases  in  strength.  The  basic 
scaling  is  encompassed  in  the  QZQ  hypothesis  which  states 

= constant  (1) 

where  Q is  the  delayed  energy  (or  waste  heat)»  Z is  a radial  coordinate 
corrected  for  the  mass  of  the  explosion  and  q is  a constant  over  a 
wide  range  of  explosions  having  the  ideal  values  of  3.5  for  strong 
shocks  and  4.0  for  weak  shocks.  If  the  prompt  energy  remaining  in  the 
blast  is  known  for  any  blast  radius  then  the  entire  blast  field  can  be 
scaled  using  Equation  (1). 

The  blast  overpressure  at  any  radial  location  can  be  determined 
once  Q(Z)  is  known.  Porzel-^  derives  an  e.xpression  for  Q as  a function 
of  overpressure  based  on  the  Rankine-Hugoniot  equations  and  a 
generalized  equation  of  state  (GOS) . He  also  gives  a binomial 
expansion  for  very  low  overpressures  and  an  empirical  relation  for 
very  high  overpressures. 

The  scaling  embodied  in  Equ.'»tion  (1)  has  been  highly  successful  in 
predicting  the  blast  field  for  spherical  e.xplosions  from  one  pound  of 
TNT  to  nuclear  explosions.  Application  of  the  QZQ  hypothesis  to  gun 
muzzle  blast  is  complicated  by  two  factors:  calculation  of  the 
equivalent  explosion  strength  or  hydrodynamic  yield  and  the  high  degree 
of  asymmetry  of  the  blast  field  compared  to  a sperical  explosion. 

B.  Equivalent  Hydrodynamic  Yield  for  a Gun  Muzzle  Blast 

A central  concept  in  IJTE  is  the  division  of  the  explosion  energy 


into  prompt  energy  which  supports  the  blast  wave  and  delayed  energy 
which  does  not.  The  prompt  energy  includes  the  energy  due  to  static 
overpressure  and  dynamic  pressure  and  also  includes  the  kinetic  energy 
due  to  material  velocity.  For  a gun  the  prompt  energy  available  to 
support  the  blast  is  a small  fraction  of  the  total  propellent  energy. 

A significant  fraction  of  the  propellent  energy  is  used  to  impart 
kinetic  energy  to  the  projectile.  Of  the  remaining  energy  a signifi- 
cant portion  is  lost  in  turbulent  boundary  layer  generation  in  the 
gun  barrel  and  impedance  of  the  flow  by  the  tube  roughness. 

In  the  previous  application  of  UTE  to  the  muzzle  blast  problem 
Hcnriksen  and  Cummings  have  analyzed  the  barrel  loss  processes  and 
have  shown,  based  on  arguments  first  advanced  by  Porzel,  that  of  the 
remaining  energy  only  one-sixth  is  prompt  energy  which  can  support  the 
blast  waves.  In  their  analysis  Henriksen  and  Cunnings  did  not  include 
the  kinetic  energy  of  the  moving  charge  in  the  determination  of  hydro- 
dynamic  yield,  but  included  the  directed  kinetic  energy  in  their  treat- 
ment of  the  non-spherical  nature  of  the  blast  field. 

In  the  present  analysis  the  kinetic  energy  of  the  moving  gas  at 
the  instant  the  projectile  leaves  the  barrel  is  included  in  the  hydro- 
dynamic  yield.  Tne  initial  prompt  energy  is 

^ 

where  M is  the  propellent  mass,  e is  the  specific  energy  of  the  pro- 
pellent, KE  is  the  projectile  kinetic  energy,  6 is  the  barrel  loss 
factor  as  given  by  the  HC  theory  and  Vj,  is  the  muzzle  velocity.  An 
initial  value  of  the  explosion  radius  is  required  in  addition  to  the 
initial  prompt  energy  in  order  to  use  the  QZQ  hypothesis  to  predict 
blast  overpressures.  For  the  gun  muzzle  blast,  the  initial  explosion 
may  be  considered  as  the  mass  of  high  pressure  gas  which  exists  the 
muzzle  immediately  after  the  projectile  has  left  the  muzzle.  Schlieren 
photographs  taken  at  these  times  show  a "banrel  shock”  system  whose 
characteristic  size  is  of  the  order  of  two  muzzle  diameters.  Hence,  in 
the  present  analysis  the  barrel  diameter  is  taken  as  the  initial  radius. 
Further,  the  QZQ  scaling  is  relatively  insensitive  to  the  initial 
radius  chosen - 

C.  The  Non-spherical  Geometry 

Experimental  measurements  of  the  blast  fields  produced  by  guns 
have  shown  that  for  a fixed  distance  from  the  observer  to  the  muzzle, 
considerably  larger  overpressures  occur  in  the  region  foreward  of  the 
muzzle  than  in  the  region  behind  the  muzzle.  This  is  analogous  to  the 
moving  charge  effect  which  has  been  studied  by  Armendt  and  Sperrazza'^. 

4 B.  F.  Anaendt  and  J.  Sperrazza,  "Air  Blast  Measurements  Around 

Moving  Explosive  Charges,  Part  III",  BRL  Memorandum  Report  No.  1019, 
V.  S.  Armg  Ballistic  Research  Laboratories  Aberdeen  Proving  Ground, 
HD,  July,  1956. 
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They  found  that  the  blast  from  a moving  charge  remained  essentially 
spherical  about  an  origin  which  moved  with  the  center  of  mass  of  the 
decelerating  charge.  The  deceleration  of  the  charge  could  be  determined 
by  conservation  of  momentum.  The  center  of  mass  decelerates  rapidly 
as  the  expanding  shock  wave  engulfs  an  ever-growing  mass  of  initially 
stationary  air. 

Detailed  calculations  were  made  based  on  this  effect  and  while  it 
predicts  higher  overpressures  forward  of  the  muzzle  than  aft,  the 
effect  is  too  small  to  explain  the  order  of  magnitude  differences  which 
are  observed  in  the  muzzle  blast  experiments.  This  is  due,  primarily, 
to  the  rapid  deceleration  of  the  center  of  mass  which  is  produced  by  a 
rapidly  expanding  spherical  shock. 

Porzel^  described  a concept  called  generalized  divergence  (GDV) 
which  permits  an  extension  of  UTE  to  non-spherical  geometries.  He 
notes  that  the  physical  significance  of  the  spatial  coordinate  in  the 
hydrodynamic  equations  is  a radius  of  curvature  rather  than  the  loca- 
tion relative  to  some  earlier  position.  It  is  the  local  radius  of 
curvature  of  the  wave  front  which  determines  its  divergence.  Thus  it 
is  the  initial  shape  of  the  charge  which  determines  the  shape  of  the 
blast  field. 

The  initial  charge  for  a gun  muzzle  blast  is  the  barrel  gases 
which  exit  when  the  projectile  leaves  the  muzzle.  The  high  over- 
pressure of  these  gases  causes  a "barrel  shock"  system  to  form  as  shown 
schematically  in  Figure  1.  This  provides  an  initial  source  for  the 
blast  field  which  is  far  from  spherical.  The  local  radius  of  curvature 
is  much  greater  in  the  foreword  portion  of  the  shock  system  resulting 
in  a smaller  divergence  and  less  rapid  decrease  in  the  blast  over- 
pressure forward  of  the  gun  muzzle  cohered  to  the  rear  where  the 
radius  of  curvature  is  much  smaller I 

The  concept  of  GOV  was  incorporated  into  the  scaling  of  blast  over- 
pressure provided  by  the  QZQ  hypothesis  by  modifying  the  distance  scale 
at  each  angular  position  in  the  blast  field  to  correct  for  the  initial 
non-spherical  gccxaetiy.  The  value  of  R used  in  the  QZQ  calculation  was 
given  by 


R = R(0)/G(8) 


(5} 


where  G(6)  is  a geometry  factor  which  was  determined  on  an  ad  hoc  basis. 
It  was  found  that  the  experimental  data  were  reasonably  well  repre- 
sented using  a geometry  factor  given  by 


G(0)  = ^Icos  0 + r^cos^^0  ♦ (4/3)^] 


(4) 


The  derivation  of  this  function  was  based  on  a velocity  argtiaent  and  is 
given  in  Appendix  A. 

To  summarize,  the  calculation  of  overpressure  for  gun  muzzle 
blast  was  accomplished  by  the  following  sequence: 


-SB, I- 


I 
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1.  Calculate  the  initial  hydrodynamic  yield  using  Equation  (2). 

2.  T!»e  hydrodynamic  yield  was  used  to  determine  the  constant  in 
the  scaling  law.  Equation  (1). 

5.  At  a given  position  (distance  and  angle)  in  the  blast  field 
the  value  of  Q was  computed  from  Equation  (1)  using  a value  of 
R determined  from  Equations  (3)  and  (4) . 

4.  The  blast  overpressure  at  the  given  field  position  was 

determined  from  the  value  of  Q determined  in  step  3 using  the 
relations  given  by  Porzel^. 

t).  Pulse  Length  Calculations 

Henriksen  and  Cummings  determined  the  pulse  length  of  the  over- 
pressure pulse  based  on  a characteristic  length  and  the  particle  velo- 
city. The  characteristic  length  was  representative  of  the  volume 
occupied  by  the  remaining  prompt  energy  in  the  blast. 


An  alternate  method  of  calculating  the  pulse  length  has  been 
developed  based  on  the  propagation  of  waves  of  finite  amplitude.  As 
the  blast  wave  from  the  e.xplosion  expands,  the  initial  high  over- 
pressure at  the  shock  front  decreases  behind  the  shock  front.  At  some 
point  in  the  c.xpansion,  designated  the  transition  radius,  the  over- 
pressure behind  the  shock  drops  to  zero  and  upon  further  expansion  the 
blast  wave  develops  a negative  phaseS.  The  fluid  velocity  behind  the 
shock  has  a similar  behavior  with  a large  velocity  in  the  direction  of 
the  shock  velocity  decreasing  to  zero  and  becoming  negative.  The 
resulting  pulse  forms  are  shown  in  Figure  2.  These  pulses  of  finite 
amplitude  propagate  into  the  undisturbed  air.  The  propagation  velocity 
is  not  the  same  for  all  portions  of  the  pulse,  but  depends  on  the  local 
speed  of  sound  and  the  fluid  velocity.  The  leading  edge  of  the  pulse 
propagates  with  a greater  velocity  than  the  rest  of  the  pulse  and 
therefore  the  pulse  increases  in  length  as  it  propagates  outward. 


The  local  wave  speed  for  the  pulse  of  finite  amplitude  is  given 

byo 


c = a^  + u (5) 

where  a and  u are  the  local  value  of  the  soecd  of  sound  and  fluid 
n * 

velocity.  For  an  isentropic  process  the  local  wave  speed  becomes 

c = a ♦ u (6) 

where  a is  the  speed  or  sound  in  the  undisturbed  region  ahead  of  the 


5.  J'u.B.  Zel’dovich  and  Yu.  P.  Raizer,  Physics  of  Shock  Waves  and  High- 
Temperature  Htidrodynamic  Phencmena,  Vol.  1 Academic  Press,  New  York, 
1966. 

6.  a.  hi.  Liepmann  and  A.  Roshko,  Elements  of  Gasdynamics,  John  Wiley 
and  Sons,  New  York,  1957. 
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pul sc ^ The  propagation  of  a shock  is  not  an  isentropic  process,  but 
for  overpressures  at  distances  greater  than  the  transition  radius  the 
corrections  are  SBiall  compared  to  the  uncertainty  in  the  measured  val- 
ues. 


The  trajectories  of  the  pulse  front  and  the  point  at  which  the 
fluid  velocity  drops  to  zero  are  shown  in  Figure  3.  The  pulse  shape  is 
shown  at  the  transition  radius  r^  and  at  two  ether  positions. 

Ne  will  define  the  pulse  length  t as  the  time  between  these  two 
trajectories  at  a given  position  r.  Note  that  as  the  shock  wave 
expands  and  the  fluid  velocity  behind  the  shock  approaches  zero,  the 
shock  front  approaches  the  speed  of  sound  and  the  pulse  length 
approaches  an  asv'nqitotic  value  which  then  propagates  as  an  acoustic 
wave. 


The  zero-velocity  trajectory  is  given  by 


^1  = 


r ) ♦ t ♦ T 
o o o 


and  the  shock  front  velocity  by 

dx 


- ” 3r  I • 2 


Y*1 

2 


— ♦ t 
u o 


(7) 

(S) 


TIjc  difference  of  these  expressions  gives  the  pulse  length 


T = T .1 

o a 


(r-to) 


r - ' - 

K'  .7 


where  r is  the  initial  pulse  length  at  the  transition  radius  r . 

I o__  ^ ^ ^9  ••  _o« 


(9) 


The 


initial  pulse  length  cai.  be  estimated  with  the  aid  of  Figure  4 which 
shows  the  overpressure  at  the  time  t when  the  shock  reaches  the  tran- 
sition radius^.  This  is  the  radius  for  which  the  negative  phase  has 
fully  developed  and  occurs  at  a pressure  ratio  of  two  across  the  shock. 
If  we  make  the  ad  hoc  assumption  that  the  pulse  occupies  half  the 
spherical  radius  at  transition  then  the  characteristic  initial  pulse 
length  is  given  by 


T = 


O 

O 2Cj; 

where  is  the  wave  speed  of  the  shock  at  a pressure  ratio  of  twc. 


(10) 


The  pulse  length  is  calculated  from  Equation  (9)  using  the  initial 
pulse  length  given  in  Equation  (10)  and  the  fluid  velocity  determined 


7.  H.  L.  Brode,  "Nuaerical  Solutions  of  Spherical  Blast  Waves",  J. 
Appl.  Phys.,  Vol.  26,  No.  6,  June  1955,  pp.  766-775. 
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frora  the  Rankine-Hugoniot  relations^  and  the  local  value  of  over- 
pressure determined  frcM  UTE. 

III.  RESULTS 

The  purpose  of  this  study  was  to  determine  the  applicability  of 
the  Hcnriksen  and  Cuanings  theory  for  smaller  caliber  guns  and  to 
extend  the  analysis  to  include  those  positions  aft  of  the  mutzle  exit 
plane.  The  theory  has  been  applied  to  predict  nuzzle  blast  over- 
pressure for  guns  ranging  from  an  8 inch  diameter  naval  gun  to  a .30 
caliber  pistol.  The  predictions  of  the  theory  have  been  collared  to 
the  data  of  Westine®  and  the  results  are  presented  in  Figures  5 
through  9. 

The  theory  is  clearly  capable  of  predicting  the  muzzle  blast  over- 
pressure for  guns  of  vastly  different  calibers.  Careful  examination  of 
the  curves  shows  that  the  theory  will  also  predict  the  angular  vari- 
ations in  the  blast  field  with  reasonable  accuracy  over  the  full  range 
of  calibers  examined.  The  largest  discrepency  between  theory  and 
experiment  occurs  at  distances  beyond  50  calibers  where  the  experimental 
data  exceeds  the  theoretical  predictions.  The  experimental  data  were 
all  obtained  with  the  guns  firing  essentially  horizontally  over  the 
ground  or  a ground  plane.  The  shock  wave  associated  with  the  muzzle  blast 
will  reflect  from  this  plane  causing  a reinforcement  to  the  primary 
shock  which  increases  the  measured  overpressure.  Indeed » in  mich  of 
the  data  the  measured  overpressure  values  increase  with  increasing  dis- 
tance from  the  muzzle  for  distances  beyond  SO  calibers.  Since  the 
theory  does  not  account  for  reflected  shocks  the  disagreement  at  large 
distances  is  not  surprising. 

Examination  of  the  angular  behavior  of  the  theory  for  distances 
less  than  50  calibers  for  the  cases  considered  reveals  that  on  the 
average  the  theoretical  predictions  tend  to  overestimate  the  over- 
pressure in  the  aft  (greater  than  90*)  portion  of  the  blast  field  and 
underestimate  the  overpressure  in  the  foreward  portion  of  the  field. 
This  is  particularly  noticeable  for  the  .30  caliber  pistol,  the  smallest 
gun  considered.  This  effect  is  a result  of  the  particular  form  of  the 
geometry  factor  chosen.  An  iaprovement  might  be  obtained  by  modifying 
the  form  of  the  geometry  factor  according  to  the  ratio- of  random  to 
ordered  motion  in  the  nuzzle  gases  as  proposed  in  Reference  2. 

Since  the  current  theory  contains  the  kinetic  energy  of  the 
mizzle  gas  in  the  prompt  energy  the  overpressure  was  calculated  for  a 
105  mm  gun  in  which  the  muzzle  velocity  varied  as  a result  of  different 
propellent  charges.  The  results  are  shown  in  Figure  10.  The  theory 

8.  P.  S.  Vestine,  "Ifodeling  the  Blast  Fielil  Around  Naval  Guns  and 
Conceptual  Design  of  a Model  Gun  Blast  Facility,"  TR  02-2643-01 
Southwest  Research  Institute,  Septeaber,1970,  AD  875  984. 
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accurately  predicts  mizzle  blast  overpressure  for  these  cases  although 
the  Mizzle  velocity  is  higher  by  a factor  of  tvo  for  the  zone  S (ZS) 
co^ared  to  the  zone  5 (Z5)  experiment.  The  result  shows  that  the 
theory  will  predict  the  Mizzle  blast  from  guns  over  a considerable 
range  of  mizzle  velocity. 

I*ulsc  length  of  the  overpressure  pulse  is  M>re  difficult  to 
CoWare  since  there  arc  large  variations  in  the  experimental  values 
obtained.  A t>'pical  situation  is  shown  in  Figure  11  for  the  case  of  thf 
same  20  am  gun  shown  in  Figure  7.  The  theory  predicts  that  the  pulse 
length  will  increase  with  distance  from  the  muzzle  until  it  reaches 
seme  as>'iq)totic  value  but  the  experiments  often  show  that  the  pulse 
length  first  increases  then  decreases  with  distance.  This  decrease 
first  appears  at  distances  greater  than  SO  calibers  and  aay  be 
related  to  the  shock  reflection  from  the  ground  plane  which  produces 
pulses  of  cc^lex  shapes. 

Predicted  pulse  lengths  for  four  guns  ranging  from  an  S inch  naval 
gun  to  a .30  caliber  rifle  arc  shown  in  Figure  12.  The  pulse  length 
was  calculated  for  an  angle  of  90*  to  the  line  of  fire  and  the  experi- 
acntal  data  were  taken  at  73*  and  103*.  The  measured  values  of  pulse 
length  are  reasonably  well  predicted  by  the  theory  over  the  wide  range 
of  pulse  lengths  prcwhiced  by  guns  of  greatly  different  caliber. 

These  results  indicate  that  the  current  theory  gives  reasonable 
agreeaent  with  experiment  over  a wide  range  of  gun  calibers.  The 
theory  has  the  advantage  of  siaqilicity  and  ease  of  calculation  and 
should  be  of  great  value  in  the  prodicticn  of  free-ficld  gun  mizzle 
blast  effects. 


l'»i|IMll|!illll|l»l|lllll|l|lllllll|l|lllll|l|||||lllll|l||l||l||l|l||lll|l||^^  |l  
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FIGURE  2. 


Sketch  of  Overpressure  and  Fluid  Veil 
After  the  Negative  Phase  lias  Developi 
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Pressure  ratio 


FIGURE  4.  Overpressure  Pulse  at  Transition.  The  Pressure  Ratio 

Drops  to  Unity  of  Approximately  One-Half  the  Transition 

Radius  r . 
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FIGURE  5.  Overpressure  from!  8 Naval  Gun 
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FIGURE  6.  Overpressure  from  3 inch  Naval  Gun. 


figure  7,  Overpressure  from  20  mm  Gun 


Overpressure  from  .30  Caliber  Rifle. 


Pistol . 


FIGURE  11.  Overpressure 


APPENDIX  A 


The  form  of  the  geometr)'  function  G{e)  was  developed  on  the 
basis  of  a moving  charge.  If  the  velocity  at  soro  point  in  the  blast 
field  is  the  sum  of  ^he  charge  velocity  V^.  and  the  blast  velocity  for 
a stationary  charge  Vg  then  V = is  assumed  to  be  in  the 

direction  of  the  radius  vector  from  the  muzzle.  If  the  ficjd  position 
is  determined  ijv  the  distance  from  the  muzzle  and  the  pol^r  angle 
relative  to  the  line-of-fire  then  6 is  the  angle  between  V and  Vj.  and 

V"  = V“  + V“  - 2VV  cose  A1 

s c c 


V V '^c 

I ■ v~  * 
si  s s 


Solving  for  the  ratio  V/V^  we  obtain 


-1=0 


^ jeosO  ± lcos-9  - [ 1 - 

h'e  assume  that  the  gcometiy  factor  has  a similar  form  ic. 
G{0)  = a ^cosG  + ^cos“e  + 6^  ^ 


G{0)  = a 


where  a and  8 arc  constants  to  be  determined. 

If  we  require  that  at  the  angle  e=ii/2  the  gecmictr)'  factor  have  the 
value  unity  then 

8 = 1/a 

licnricksen  and  Cummings'  find  that  the  kinetic  energy  of  ordered 
motion  at  the  muzzle  is  approximately  twice  that  of  random  motion. 
Based  on  this  result  we  require 

G(0)  = 2 


This  gives  the  results 
a = 3/4 


G(e)  = j jocose  + f cos"0 


”7] 
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ABPENDJX  B r 

A BASIC  coBpiitJzr  *progi:^  Ka4  Kcim  cpnsitructud  to  pcrfora  the 
calculations  for  tlie- euzzIc  blast  field  based  on  the  preceding  analysis. 
The  prograa  is  listed  at  the  end  of  this  appendix,  llic  input  para^tcis 
required  for  the  prograa  are  as  folloKs: 

1.  Projectile  aass  - kilograas  (kg) 

2-  lizzie  velocity  - actres/second  (a/s) 

3.  Propellent  specific  energy  - joules/kilograa  (.l/kg) 

4.  Ratio  of  peak  chaeber  pressure  to  ataospheric  pressure  - 
diacnsionless 

5.  Barrel  length  - setres  (a) 

6.  Barrel  diaacter  - metres  (a) 

7.  Height  of  the  barrel  rifling  - aetres  (a) 

hhen  execution  of  the  program  begins  the  program  will  pause  for  input 
until  these  values  are  entered  in  the  order  given  above,  hben  the  last 
value  is  entered  the  prograa  calculates  the  overpressure  at  the  choke 
position  and  at  the  muzzle  and  an  equivalent  explosive  yield.  Next, the 
program  will  pause  to  allow  optional  values  of  Ql,  Q2  and  ratio  of 
specific  heats,  y to  be  entered  if  desired.  The  variables  Q1  and  Q2 
arc  the  exponents  in  the  Porzcl  QZQ  h>*pothcsis  for  the  strong  and  weak 
shock  regimes  and  have  the  standard  values  3.25  and  5^5  respectively. 
The  standard  value  of  y is  1.4.  .After  the  lotion  is  exercised  the 
program  co^utes  the  transition  radius  and  pauses  until  the  desired 
angular  position  in  the  blast  field  is  entered.  This  is  the  angle  @ 
measured  from  the  line  of  fire  to  the  radius  vector  fr<m  the  aizzlc  to 
the  field  location  and  lies  between  0*  and  ISO*. 

The  program  is  designed  to  calculate  overpressure  and  pulse 
length  as  a function  of  distance  fr«»:  the  muzzle  at  the  input  angle 
0.  The  prograa  will  pause  until  the  desired  distance  closest  to  the 
EUizzle  is  entered.  It  then  pauses  until  an  increment  in  the  radial 
distance  is  entered  and  again  inttil  the  total  number  of  increments 
desired  is  entered.  The  prograa  then  cfsputcs  the  overpressure  and 
pulse  lengths  at  each  radial  location  specified  for  the  angle  S 
specified,  hhen  these  co^utations  are  co^leted  the  prograa  pauses 
to  allow  the  option  of  terminating  the  prograa  or  specifying  another 
angle  8^  If  a new  angle  is  entered  then  the  prograa  pauses  to  allow 
an  option  of  continuing  with  the  saae  radial  locations  or  changing 
to  new  radial  locations.  Khcn  the  last  angular  position  desired  has 
been  calculated  the  operator  enters  ainus  one  (-1)  to  terminate 
execution. 
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